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Abstract. - We investigate the dynamics of polymer translocation through nanopores under 
external driving by 3D Langevin Dynamics simulations, focusing on the scaling of the average 
translocation time r versus the length of the polymer, r ~ N a . For slow translocation, i.e., under 
low driving force and/or high friction, we find a ~ 1 + v ~ 1.588 where v denotes the Flory 
exponent. In contrast, a ~ 1.37 is observed for fast translocation due to the highly deformed 
chain conformation on the trans side, reflecting a pronounced non-equilibrium situation. The 
dependence of the translocation time on the driving force is given by r ~ F" 1 and r ~ p-oso £ Qr 
slow and fast translocation, respectively. These results clarify the controversy on the magnitude 
of the scaling exponent a for driven translocation. 
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Introduction. — The passage of a polymer through 
a nanopore is essential for numerous biological processes 
such as DNA and RNA translocation across nuclear pores, 
protein transport through membrane channels, or virus in- 
jection into cells [1]. (Bio)polymer translocation is also at 
the heart of various potential technology applications, for 
instance, rapid DNA sequencing, gene therapy, or con- 
trolled drug delivery [2]. A translocating polymer has to 
overcome an entropic barrier. In biological cells, biopoly- 
mcrs are driven through pores by transmembrane poten- 
tials, chemical potential gradients due to binding proteins, 
active pulling by polymerase, or entropic pressure (virus 
ejection) [1]. From both the basic physics as well as a 
technology design perspective, an important measure is 
the scaling of the average translocation time r with the 
polymer length N, t ~ N a , and the value of the corre- 
sponding scaling exponent a. 

Most translocation experiments are carried out under 
an applied electric field across the pore. We focus here 
on this particular type of driven translocation. In exper- 
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iments with a-hemolysin pores of inner diameter 2 nm a 
linear behavior of r ~ N (a = 1) was observed [3,4], while 
an exponent a = 1.27 was obtained for double-stranded 
DNA translocation through a solid-state nanopore of in- 
ner diameter 10 nm [5]. Recently the voltage-driven 
translocation of individual DNA molecules through solid- 
state nanopores of diameters 2.7 ... 5 nm revealed that 
t ~ N 1A0 for DNA molecules in the range 150 - 3500 
base pairs [6]. 

Inspired by these experiments, a number of recent the- 
ories and simulations on the translocation dynamics have 
been presented. Standard Kramers analysis of diffusion 
across an equilibrium entropic barrier yields r ~ N 2 for 
unbiased translocation and r ~ N under external driving 
(assuming friction to be A-independcnt) [7,8]. However, 
as noted in Ref. [9] the quadratic scaling behavior in the 
unbiased case cannot be correct for a self-avoiding poly- 
mer, as the translocation time would be shorter than the 
Rouse equilibration time scaling like tr ~ N 1+2v , involv- 
ing the Flory exponent v (y = 0.588 in 3D, v = 0.75 
in 2D) [10]. This finding renders the concept of equilib- 
rium entropy and the ensuing picture of entropic barrier 
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crossing inappropriate for translocation dynamics. Nu- 
merically, it was shown that for large N, t ~ N l+2u , i.e., 
the translocation time scales in the same way as the equi- 
libration time, but with a much larger prefactor [9]. This 
result was recently corroborated by extensive numerical 
simulations based on the Fluctuating Bond (FB) [11] and 
Langcvin Dynamics (LD) models with the bead-spring ap- 
proach [12-14]. For driven translocation, the deviation 
from equilibrium is expected to be even more pronounced, 
and the equilibrium cntropic barrier is less relevant for 
the translocation dynamics. For this lower bound 

t ~ N 1+l/ was estimated for the translocation time [15]. 
Lattice Monte Carlo (MC) simulations of self-avoiding 
chains in 2D revealed a « 1.5 [15], which is smaller than 
1 + v = 1.75, a difference attributed to finite size effects. 
However, additional simulation studies found a crossover 
from r ~ jyi-46io.cn w jy2v f or re i a y V ely short polymers 
to r - jyi-70±o.03 w ^1+^ f or i on gcr chains (N > 200) us- 
ing FB [16] and LD [12, 17] models in 2D. Increasing the 
friction, the crossover vanished and only exponent 1 + v 
was observed [12]. However, the crossover is absent in 3D 
for quite long chains (N ~ 800) and a rs 1.40 is observed 
using both LD [18] and GROMACS molecular dynamics 
with LD thermostats [19] simulations. Lchtola et al. [20] 
find approximately the same exponent of about 1.40 based 
on LD simulations, however, their exponent increases with 
increasing force, which cannot hold asymptotically. More- 
over, Gauthier and Slater [21] reported r ~ N 1+L/ using 
an exact numerical method valid at low bias. 

Recently, however, alternative scaling scenarios have 
been presented in Refs. [22, 23] that contradict above re- 
sults. These two views disagree with each other [24]. To 
resolve the apparent discrepancy on the value of a for 
driven translocation, we here report results on the be- 
havior of t as a function of N from high-accuracy LD 
simulations in 3D. As in the above-mentioned studies, we 
also neglect hydrodynamics effects as well as polymer-pore 
interactions [25]. We find that there exists two limit- 
ing regimes, corresponding to slow and fast translocation 
respectively. The slow translocation case is realized for 
low driving forces and/or high friction, and in this regime 
a « 1 + v. In the opposite limit of fast translocation for 
high driving forces and/or low friction, the corresponding 
scaling exponent is given by a « 1.37. As we will argue 
below, the difference between the scaling exponents for 
these two regimes can be ascribed to the highly deformed 
chain conformation during fast translocation, reflecting a 
pronounced non-equilibrium situation. Our results clarify 
the controversy on the value of a for driven translocation. 

Model and method. — In our simulations, the poly- 
mer chains are modeled as bead-spring chains of Lcnnard- 
Joncs (LJ) particles with the Finite Extension Nonlinear 
Elastic (FENE) potential. Excluded volume interaction 
between monomers is modeled by a short range repulsive 
LJ potential: U LJ (r) = 4e[(f ) 12 - (f ) 6 ] +e for r < 2^ 6 a, 
and for r > 2 1 / 6 ct. Here, a is the monomer diameter 




Fig. 1: Schematic representation of our simulation system. The 
wall consists in a single layer of beads while the pore itself is 
formed by eight particles with their centers equally distributed 
on a circle of diameter 3a. Therefore, the actual pore diameter 
is 2a. 

and £ is the potential depth. The connectivity between 
neighboring monomers is modeled as a FENE spring with 
Ufen e(t) = -^fci?o m (l — r2 /^o)i where r is the dis- 
tance between consecutive monomers, k the spring con- 
stant, and R the maximum allowed separation between 
connected monomers. The wall ("membrane") carrying 
the pore is composed of particles of diameter a. The wall 
thickness is a. The nanoporc consists of eight particles 
with their centers equally distributed on a circle of diam- 
eter 3a, see Fig. [TJ Therefore, the actual pore diameter is 
2cr. Between all monomer-wall particle pairs, there exists 
the same short range repulsive LJ interaction as described 
above. 

In LD simulations each monomer is subjected to conser- 
vative, frictional, and random forces: mrj = — V(C/j,j + 
Ufene) — +F cx t +Ff [26], where m is the monomer 
mass, £ the friction coefficient, and ~Ff the random force, 
that satisfies the fluctuation-dissipation theorem. The ex- 
ternal force is expressed as F cxt = Fx, where F is the ex- 
ternal force strength exerted on the monomers in the pore, 
and x is a unit vector in the direction along the pore axis. 
In the present work, we use the LJ parameters e and a, and 
the monomer mass m to fix the energy, length and mass 
scales. This sets the time scale ti,j = (ma 2 /e) 1 ^ 2 . The 
dimensionless parameters in our simulations are Rq = 2, 
k = 7, fcsT = 1.2 and £ = 0.7, unless otherwise stated. 
The Langevin equation is integrated in 3D by a method 
described in Ref. [27]. Initially, the first monomer of the 
chain is placed in the entrance of the pore, while the re- 
maining monomers evolve in the Langevin thermostat to 
obtain an equilibrium configuration. The translocation 
time is defined as the time interval between the entrance 
of the first segment into the pore and the exit of the last 
segment. Typically, we average our data over 1000 inde- 
pendent runs. 

We note that in our model hydrodynamic interactions 
and the polymer-pore interactions are neglected. Regard- 
ing the issue of hydrodynamics, recent Molecular Dynam- 
ics [28] and Lattice Boltzmann [29,30] results show that 
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Fig. 2: Translocation time r versus chain length N for different 
driving forces F and friction coefficients £. For fast transloca- 
tion (low driving force or high friction) a w 1.37; for slow 
translocation tissl + i/, 



hydrodynamics is screened out in a narrow pore, which 
is the case studied here as well as in the experiments. 
In Ref. [31], however, minor increases of scaling expo- 
nents with increasing driving force are obtained in the 
presence of hydrodynamic interactions. Recently, we have 
used LD simulations to investigate the influence of an at- 
tractive component of polymer-pore interactions on the 
translocation dynamics [25]. We found that with increas- 
ing strength of the attractive interaction, the histogram 
for the translocation time r shows a transition from a 
Gaussian distribution to a long-tailed distribution corre- 
sponding to thermal activation over a free energy barrier. 
In addition, a strong attractive polymer-pore interaction 
can directly affect the effective scaling exponents of r both 
with N and with the applied voltage, which provides a pos- 
sible explanation for the different experimental findings on 
these physical quantities. However, our main focus here is 
to arrive at a detailed understanding of the translocation 
dynamics for a purely repulsive pore, and to compare with 
the results in Refs. [22, 23] to settle the controversies re- 
garding the scaling behavior of translocation time versus 
the length of the polymer. Hence, the attractive part of 
the polymer-pore interaction is left out in this study. 

Results and discussion. — As defined above, our 
pore length is a and the pore diameter is 2a. Previous 
studies [16] showed that for longer pores there is no ob- 
vious power-law scaling for relatively short chains. This 
may be the reason for a linear dependence reported in 
Ref. [32] where the pore length is 5er with chain length 
N < 80 and in Ref. [33] where the pore length is 12cr 
with chain length N < 100. Various heuristic scaling ar- 
guments for r have been presented for short pores, e.g., 
in Refs. [15, 16, 22, 23, 34] and will be compared to our 
numerical results below. 

The translocation time as a function of the polymer 
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Fig. 3: Scaling of translocation velocity with chain length. 



length is plotted in Fig. [2] One of the main features is 
that a depends significantly on both driving force and 
friction. For stronger driving forces (F = 10.0 and 5.0) 
and lower friction (£ = 0.7), we find a = 1.37 ± 0.02 and 
a = 1.37 ± 0.05, which are in agreement with previous re- 
sults [18,19,35] and the recent experimental findings for in- 
dividual short DNA molecules in the range 150 ~ 3500 bps 
through solid-state nanopores [6]. However, decreasing F 
or increasing £, we observe a = 1.58 ±0.03, 1.51 ±0.01 
and 1.52 ±0.01 for F = 0.5 and £ = 0.7, F = 2.5 and 
£ = 3.0, and F = 5.0 and £ = 3.0, respectively. These 
exponents arc in good agreement, not only with the pre- 
diction l + i/ from Ref. [15,34], but also with the results 
from the exact numerical method reported in Ref. [21] for 
low fields, as well as with our previous 2D simulations 
for relatively long polymers [12,16,17]. These observa- 
tions demonstrate that there exist two limiting dynamic 
regimes as demonstrated in Fig. slow and fast translo- 
cation. For slow translocation, a ks 1 + v, while a ~ 1.37 
for fast translocation. 

In the scaling arguments presented in Ref. [15] , an essen- 
tial assumption is that the chain is not severely deformed 
during translocation, and in particular the chain configu- 
ration of the already translocated part of the chain is close 
to equilibrium. For slow translocation, these assumptions 
are satisfied. Thus, it is not surprising that the predicted 
exponent 1 + v is observed in this regime, in contrast to 
the fast dynamics regime claimed in Ref. [20] . Conversely, 
for fast translocation the chain is highly deformed and the 
translocation dynamics is different. In our previous 2D 
simulations based on FB [16] and LD [12,17], a crossover 
from t ~ N 2v for relatively short polymers (N < 200) to 
t ~ N 1+v for longer chains was found under a stronger 
driving force F = 5 and friction £ = 0.7. Increasing N also 
slows down the translocation dynamics and thus the sys- 
tem changes from the fast to the slow translocation regime. 
However, for the same F and £ values we failed to observe 
a similar crossover for N = 40 ~ 800 in 3D [19], possibly 
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due to the fact that the value of N for the crossover is at 
a much higher value as compared to the 2D situation. To 
actually observe the crossover, an alternative choice is to 
lower the driving force or increase the friction. 

Based on the fractional Fokker-Planck equation involv- 
ing long-range memory effects (compare Ref. [36]), the 
scaling a = 2v + 1 — 71 was obtained, such that a = 1.55 
in 2D and 1.5 in 3D; here 71 (« 0.945 2D and w 0.68 in 
3D) is the critical surface exponent [23]. The MC simula- 
tions in Ref. [23] were carried out for weak driving forces, 
producing r ~ N /F in 3D, in excellent agreement with 
the value t ~ N l+L> /F for slow translocation processes. 
Moreover, using linear response theory with memory ef- 
fects [24], Vocks et. al. [22] came up with an alternative 
estimate r ~ A^t^ 77 " for 3D, which means a = 1.37 in 3D. 
Their a in 3D is consistent with our numerical data for fast 
translocation. However, their estimate fails to capture the 
scaling exponent for slow translocation. 

For the driven translocation experiments with a voltage 
applied across the pore, the force F acts only on the few 
monomers inside the pore. As a consequence, it was ar- 
gued [15] that the center of mass of polymer should move 
with a velocity v ~ F/N . Thus the lower bound for the 
translocation time of an unhindered polymer is the time 
to move through a distance R g (radius of gyration of the 
polymer), giving rise to the scaling r ~ R g /v ~ N 1+v /F. 
Fig. [3] shows the translocation velocity v as function of 
the polymer length for different driving forces and friction 
coefficients, v ~ N s . Here v is the average velocity with 
respect to the last monomer (for details see Ref. [16]), and 
we checked that the corresponding v of the center of mass 
scales in the same way. For slow translocation, wc find 
5 = -1.01 ±0.02 for F = 0.5 and £ = 0.7, 5 = -0.95±0.02 
for F = 2.5 and £ = 3.0, and S = -0.94 ±0.01 for F = 5.0 
and £ = 3.0, which are in good agreement with v ~ N . 
For fast translocation the velocity decreases less rapidly 
with N: 5 = -0.79 ± 0.01 for F = 10.0 and £ = 0.7 and 
5 = -0.79 ± 0.02 for F = 5.0 and £ = 0.7, as observed 
in Ref. [18]. Based on the values of F/£ we roughly dis- 
tinguish the slow and fast regimes, see Table 1. For fast 
translocation with £ = 0.7, F/£ = 14.28 for F = 10.0 and 
F/£ = 7.14 for F = 5.0, which are much higher than those 
for the slow translocation where F/£ = 0.71 for F = 0.5 
and £ = 0.7, F/i = 0.83 for F = 2.5 and £ = 3.0, and 
F/£ = 1.67 for F = 5.0 and £ = 3.0. 

Following the change of the scaling exponent a, the 
scaling of the average velocity changes from v ~ _/v -0 - 79 
for fast translocation, to v ~ A^ 1 in the case of slow 
translocation. Furthermore, we have also studied how the 
translocation time varies as a function of the driving force 
for £ = 0.7 and N = 128 (Fig. gj. As long as F < 2 we 
observe r ~ 1/F, before crossing over to t ~ ^-0.80 f or 
strong driving. The data can be fitted with an empiri- 
cal function r = ^OF" 1 ^ + (F/2.5) 8 ] 1 / 40 , yielding the 
details of the crossover from the slow translocation (weak 
force) to the fast translocation (strong force) regime for 
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Fig. 4: Translocation time r versus driving force at £ = 0.7. 
The solid line shows the empirical fitting function for the data, 
as discussed in the text. 

this case. The seemingly high exponent 8 in the empiri- 
cal function is necessary to account for the relatively fast 
turnover between the scaling r ~ F^ 1 at low force and 
the behavior r ~ F~°- 8 at high force. The difference of 
the translocation dynamics in the two regimes can be un- 
derstood by inspecting the polymer configurations during 
the translocation process. For faster translocation, only 
part of the chain on the cis side can respond immedi- 
ately, while the remaining part near the chain end does 
not feel the force yet. As a result, a part of the chain on 
the cis side is deformed to a trumpet and even stem-and- 
flower shaped [35], while the translocated portion on the 
trans side has a compact spherical shape, as it does not 
have time to diffuse away from the pore exit, see Fig. [5] 
With increasing time, the tension propagates along the 
chain, which changes the chain conformation progressively. 
Thus, the chain cannot achieve a steady state with average 
velocity v ~ A^ 1 even in the late stage of the transloca- 
tion process. Fig. [HI shows the radius of gyration (R g ) at 
the moment just after translocation. For fast transloca- 
tion, this obeys the scaling behavior R g ~ A" - 51 , signifi- 
cantly different from the equilibrium scaling of R g ~ Af ' 6 , 
indicating a non-equilibrium compactification of the chain 
immediately after translocation. On the other hand, the 
corresponding scaling behavior for the chain after slow 
translocation is approximately the same as the chain in 
equilibrium. 

We note that in a recent theoretical study the effect 
of a trumpet shape of the chain on the cis side, on the 
translocation dynamics was found to cause a breakdown 
of the r ~ 1/F scaling [35]. However, this theory neglects 
effects due to the compacted chain structure on the trans 
side. 

We also checked the translocation coordinate (s(t)) ~ t 13 
for N = 64 for different F and £, see Fig. 7 and Table 1. 
For fast translocation with £ = 0.7, we find (3 = 0.84 and 
0.85 for F = 10.0 and 5.0, respectively. However, for slow 
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Table 1: Summary of numerical results. Here, F is the driving force, £ the friction coefficient, a the scaling exponent of 
translocation time r as a function of the chain length N, 5 the scaling exponent of translocation velocity as a function of N, and 
j3 the scaling exponent of the translocation coordinate s as a function of time. These results clearly demonstrate two regimes. 
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Fig. 5: Typical chain conformation during fast translocation 
for A^ = 128, £ = 0.7, and F — 5.0. 3D conformations are 
projected onto the XY plane. Black: early stage. Red: Later 
stage. 
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Fig. 6: The radius of gyration (R g ) before translocation and 
at the moment just after the fast translocation (£ = 0.7 and 
F = 5.0), and the slow translocation (£ = 0.7 and F — 0.5). 



translocation = 0.71 for F = 5.0 and £ = 3.0, = 0.69 
for F = 2.5 and £ = 3.0, and = 0.64 for F = 0.5 and 
£ = 0.7. Different values for fast and slow translocation 
processes demonstrate the existence of different dynamic 
regimes. The definition (s(t)) ~ t implies that N ~ t 13 . 
Compared with r ~ N a , one obtains a0 - 1. For F = 0.5 
and £ = 0.7, a.0 = 1.01 « 1, which indicates that non- 
equilibrium effect does not seem to be severe. However, 
ol0 = 1.16 for F = 5.0 and £ = 0.7 indicating a breakdown 
of "simple" scaling due to highly non-equilibrium effect. 

Conclusions. — We have investigated the dynamics 
of driven polymer translocation through nanopores by 3D 
Langevin dynamics simulations, focusing on the scaling 
of the average translocation time r as a function of the 
polymer length N. For slow translocation, i.e., under low 
driving forces and/or high friction, we find r ~ N a with 
a ~ 1 + v. In the opposite case, we obtain a w 1.37. As a 
function of the driving force F, the dependence r ~ \/F 
and t ~ ^-0.80 are obtained respectively for slow and fast 
translocation. The different behavior in the two regimes 
can be understood from analysis of the chain conforma- 
tions during the translocation process. In the slow translo- 
cation case, the configurations at all times are close to the 



equilibrium case while for the fast translocation regime, 
there exist highly deformed, unrclaxed chain conforma- 
tions throughout the translocation process. These results 
clarify the controversy on the value of a for driven translo- 
cation in the existing literature. 
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